Microscopic origin of large negative magneto-electric coupling in Sri/ 2 Bai/ 2 Mn03 
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With a combined ab initio density functional and model Hamiltonian approach we establish that 
in the recently discovered multiferroic phase of the manganite Sr^Ba^MnOs the polar distortion 
of Mn and O ions is stabilized via enhanced in-plane Mn-O hybridizations. The magnetic superex- 
change interaction is very sensitive to the polar bond-bending distortion, and we find that this 
dependence directly causes a strong magnetoelectric coupling. This novel mechanism for multifer- 
roicity is consistent with the experimentally observed reduced ferroelectric polarization upon the 
onset of magnetic ordering. 
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Multiferroic materials are ideal candidates for the re- 
alization and practical use of strong magnetoelectric ef- 
fects [lj,|2|. The scarcity of actual materials that are mag- 
netic ferroelectrics appears to be related to the compe- 
tition between the conventional mechanism of ferroelec- 
tric cation off-centering, which requires empty d-orbitals, 
and the formation of magnetic moments which requires 
partially filled d-orbitals Q • A concomitance of mag- 
netism and ferroelectricity then has to rely on more sub- 
tle microscopic coupling mechanisms, driven by spin- 
orbit coupling in the form of Dzyaloshinskii-Moriya in- 
teractions Q or exchange-striction Q • The recently syn- 
thesized manganite Sr^Bai^MnOa however defeats the 
generic incompatibility of a cation both having a mag- 
netic moment and being ferroelectrically displaced. This 
system is a classic example of a material in which charge, 
spin, lattice and orbital degrees of freedom are strongly 
coupled, giving in this particular case rise to a strong 
magnetoelectric (ME) effect, the origin of which we set 
out to clarify here. 

For doing so, the methods from modern ab initio band- 
structure theory are powerfull tools - very helpful not 
only in predicting new multiferroic materials, but also in 
understanding the underlying mechanisms for magneto- 
electric couplings. The computed values of macroscopic 
polarization P agree exceptionally well with those ob- 
served experimentally |5|— }l0| ■ In the last few years several 
ab initio calculations have pointed out the possible ferro- 
electric state with large polarization for AMnC>3, where 
A is an alkaline earth element. The proposed mecha- 
nism is based on off-centering of Mn 4+ ions stabilized 
via a charge-lattice coupling of Peierls type [lll - ll3fl . The 
problems in synthesizing such a material with predicted 
ferroelectricity has very recently been overcome: last year 
Sr 1 / 2 Ba 1 / 2 Mn03 (SBMO) has been reported to support a 
ferroelectric phase via the off-centering of magnetic Mn 4+ 
ion in conjunction with a perovskite tetragonal struc- 



ture [141] . The onset of the low-temperature long-range 
antifcrromagnetic (AFM) ordering strongly reduces the 
polarization indicating a large magnetoelectric effect [3] ■ 
The AFM order, in other words, does not support ferro- 
electricity, but it neither completely destroys it. This spe- 
cial feature of SBMO opens a new avenue for the quest of 
materials with strong ME effects, where the search need 
not be restricted to systems in which FE and magnetism 
mutually stabilize each other. 

Here we establish with a combination of first-principles 
calculations and a model Hamiltonian analysis that the 
ferroelectric polarization mainly arises from a polar dis- 
tortion of Mn and O ions caused by an enhanced in-plane 
Mn-0 hybridization. Since the magnetic superexchange 
interaction strongly depends on this distortion, a strong 
and novel type of magnetoelectric coupling arises. This 
ME coupling is negative in the sense that the ferroelectric 
polarization is not promoted by magnetism, but rather 
reduced by it, which renders antiferromagnetic ordering 
and ferroelectricity strongly coupled. 

We first present the results of our first-principles cal- 
culations based on density functional theory (DFT) 15 
using the generalized gradient approximation (PBE) 16 
and including correlation effects within the DFT+U 
scheme [l7( as implemented in VASP[l8j. We use on- 
site Coulomb and exchange parameters U=3.0 and 4.5 
eV and Jh =1-0 eV on the manganese d-orbitals. In the 
projector augmented wave scheme [l9| the cut-off for the 
plane- wave basis set was chosen as 400 eV and a 8x8x8 
mesh was used for the Brillouin-zone sampling. To cal- 
culate the electronic contribution to the spontaneous po- 
larization we use the Berry-phase method developed by 
King-Smith and Vanderbilt 2o| . In the calculations the 
in-plane lattice constant is taken as a =3. 85 A. For the 111- 
terplane distance we consider the two values c/a =1.005, 
1.01, which are experimentally determined for the SBMO 
at different temperatures 14 1. 
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FIG. 1: (Color online) Schematic view of SBMO unit cell and 
displacements of Mn and O sites in AFM magnetic struc- 
ture. Arrows indicate the relative atomic ferroelectric dis- 
placements. 



a) 



178 
177.5 

177 
176.5 

176 
175.5 



b) 




20 
10 


▲ 


P 

ionic 




-10 


■ 


Rle 


-20 
-30 


• 


P 


-40 



\ 



\ 





1= 
n 

n 



1.005 1.01 

c/a 



1.005 1.01 

c/a 



FIG. 2: (Color online) a),b) Values of a and spontaneous 
polarization P, P e i e , Pionic as function of the ratio c/a for 
17=3.0 ev and 4.5 eV respectively. 



For all the above parameters SBMO is safely in an 
AFM (G-type) insulating state with a band gap of ~ 0.4 
eV with Mn magnetic moments M ~ 2.6 fib, in agree- 
ment with previous DFT calculations on CaMnOa and 
SrMn0 3 [II [H, H3 and experiments on SrMn0 3 [Hj]. 
The valence band is predominantly majority-spin Mn t2 S 
and O 2p characther with strong p—d hybridization while 
the conduction band is formed by Mn e g orbital and 
empty minority t2 9 states, which is consistent with Mn 4+ 
in octahedral crystal field. To find the energetically most 
stable configuration we relax the ions performing struc- 
tural optimization in a 40 atoms 2x2x2 unit cell. We 
start the relaxation from a checkerboard arrangement of 
the Sr, Ba ions and check that our results do not depend 
on this assumption. 

Even if the initial ionic structure belongs to space 
group I4/mmm (No. 139) which is centrosymmetric, for 
all our parameters the relaxed structure belongs to the 
space group I4mm (No. 107) and it breaks inversion sym- 
metry. The polar ionic displacements associated with the 
reduced symmetry are shown schematically in Fig. [1] 
The O-Mn-0 angle a (see Fig. [TJ which is 180° for ideal 
centrosymmetric structure with c/a = 1, is reduced in 
agreement with the experimentally determined low sym- 
metry structure 14j. The deviation of a from 180° as 
function of c/a is shown in Fig. [5J We now analyze 
the effect of this reduced angle on the ferroelectric polar- 
ization. We first notice that the electronic contribution 
to the spontaneous polarization P evaluated in the cen- 
trosymmetric structures (I4/mmm) for the AFM ground 
state is zero, meaning that the polar state is not magnet- 
ically driven. Indeed at different values of the ratio c/a it 
is the covalent bond formation upon ionic displacements 
between e g orbitals of Mn and p orbitals of apical O ions 
to determine the stabilization of the ferroelectric state change in the superexchange interactions betweens Mn 



On the other hand in the relaxed state, the large po- 
lar displacements of the apical O ions along the c lattice 
direction result in the formation of dipolar pairs between 
manganese and oxygen (see Fig. [IJ and to a state simi- 
lar to a bond-centered charge density wave [23[. At low 
temperature, when the system orders antifcrromagneti- 
cally, this enables the practical realization of a peculiar 
and atypical multiferroic state. This is shown by the re- 
sults for the polarization P, whose electronic and ionic 
contributions (P e j e , Pionic) are plotted in Fig. [2] In- 
creasing the ratio c/a the magnitude of the electronic 
contribution P e ; e increases and that of the ionic contri- 
bution Pionic decreases. Since the two are opposite in 
sign the total polarization P = P e i e + Pionic increases as 
a function of c/a. The calculated value of P agrees with 
the experimental value of 13.5 /xC/cm 2 for single domain 
14] . This physical result does not depend strongly on the 
structural and interaction (U, Jh) parameters, but the 
quantitative description of the ferroelectric instability in 
SBMO should of course depend on the actual values of 
these parameters, as has been found to be the case in 
other Mn based multiferroic materials (f|. In particular 
at larger U the magnetic moment increases and the fer- 
roelectric tendency decreases as the angle a gets closer 
to 180° (see Fig. |2J. 

The calculated polarizations show that the FE order 
is not driven by the magnetic order, but yet the two are 
strongly coupled. This counterintuitive situation arises, 
as we will show next, from the ferroelectric transition 
in SBMO being driven by Mn and O displacement and 
the mechanism of the suppression of tetragonal distor- 
tion below Tjv [lJ] being due to the subsequent strong 
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FIG. 3: (Color online) a) Schematic view of non-collinear 
magnetic structures having Mn spins with angle 6 — 
0°,45°,90°; b) Values of a and spontaneous polarization p 
as function of the ratio c/a for different values of 9 at U=3.0 
eV and J H =1.0 eV. 



spins [24J. The 180° O-Mn-0 bonds are energetically 
favored by the antiferromagnetic coupling [24j then in 
the ferroelectric state the off-centering of Mn ions, which 
is in favor of the inset of double exchange interactions, 
gets suppressed with a net decreasing of the ferroelectric 
polarization 25] . The effect of the magnetism on the 
ferroelectric distortions can be captured by performing 
calculations with non-collinear magnetic structures hav- 
ing Mn spins with angle 9 ranging from 0° (G-type) to 
90° (see Fig. to control how the superexchange in- 
teractions along the O-Mn-0 bonds changes P. At each 
angle 9 the lattice structure is relaxed and the sum of 
electronic and ionic contributions to the ferrolectric po- 
larization is evaluated (see Fig. [3Jd). Increasing the angle 
9 between the spins reduces the superexchange intera- 
tions. The Mn-O-Mn angle a decreases with a resulting 
larger Mn off-centering which stabilizes the ferroelectric 
polarization. We observe that the magnetic order alters 
both electronic and ionic contributions to the polariza- 
tion via a change in a: the magnetism is thus coupled 
to the lattice and the latter is coupled to the polariza- 
tion. Increasing the superexchange interactions causes 
the magnetic structure to drive the lattice towards a re- 
covery of the a centrosymmetric arrangement. 

To flesh out the microscopic origin of the ferroelectric 
instability we set up a model Hamiltonian, based on the 
bandstructure results and taking into account the differ- 
ent magnetic exchange interactions and the coupling of 
the electrons to the lattice: 
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Here, d ia [pj a ) and d\ a (p\ a ) are the annihilation and 
creation operators for Mn-cZ (O-p 7 ) electrons with spin 
a =t,4- Si are the localized t^g spins (S — 3/2), which 
in this study are treated classically and coupled antifer- 
romagnetically via J s . m are the off-centering distortions 
of Mn ions along the c-axis and K s denotes the stiffness 
energy associated with these distortions. t 7 (ui) denote 
the distortion-dependent hopping amplitudes between 
rf 3z 2_ r 2 and p\ orbitals along 7 direction (7 = x,y,z). 
Note that 7 denotes the direction in real space and not 
the character of the p orbitals. The (Ji denote the elec- 
tronic spin operator defined as, erf = J2aa' 4o /7 1t<t' ^i" 1 ' > 
where r M are the Pauli matrices. A p d is the on-site en- 
ergy difference between Mn-d3 2 2_ r 2 and 0-p z levels. 

In the model Hamiltonian the ionic displacements are 
restricted to the c-axis direction, as observed in the ex- 
periments and verified in our bandstructure calculations. 
In principle the O ions are easier to displace, however a 
combination of O displacements and Mn displacements 
can be modeled as a net off-centering displacement of the 
Mn ions along with an overall change in the lattice c pa- 
rameter. Here, we model the effective displacements via 
the off-centering u^a of the Mn ions, where a is the Mn- 
Mn lattice spacing. We consider only d 3z 2_ r 2 orbital as 
the one that can hybridize with the 0-p z levels, since the 
planar orbitals d x 2_ y 2 have zero overlap with the 0-p z . 

If Mn ions are located at the center of 0$ octahedra 
then the hopping between d 3z 2_ r 2 and 0-p z is non-zero 
only along z-axis and is given by to = (pda). How- 
ever, if Mn ions are off-centered by a small displace- 
ment they lead to a finite in-plane hopping which can 
be calculated from the Slater-Koster tables as t x ' v = 
n(n 2 — {I 2 + m 2 )/2)(pda) + \/3n(l 2 + m 2 )(pdTr) where 
l,m,n are the direction cosines from O to Mn [26|. Tak- 
ing only the pda contribution one can write the hop- 
ping integral in terms of the Mn-O-Mn angle a as, 
t x /v = sin(a/2)(sin 2 (a/2) - cos 2 (a/2)/2)(pdcr). Rewrit- 
ing the trigonometric functions in terms of the the dis- 
tortions, we get to leading order in the distortion uo, 
t x /y ^ —2uq (pda). The next order term is 0(uq) which 
can be safely ignored. Naturally t z pd is also modified via 
a Peierls type term with the hopping between longer and 
shorter bonds given by t± — (1 ± gu)to. 

Just as the off-centering of Mn ions affects the hop- 
ping parameters t^ d , it also affects the value of J s via 
the Mn-O-Mn bond angle. J s is maximum at a = tt 
and is reduced by any deviation. The leading order 
change in the Taylor expansion around the point a = tt 
is 0(5a 2 ). Therefore, for small deviations we can model 
the distortion-dependence as Js^ y = — Jo cos(a). In prin- 
ciple J I is also affected since the distances Mn-Oi and 
Mn-02 for the two apical oxygens become unequal, but 
this dependence does not affect the physical picture. 
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FIG. 4: (Color online) a) u m in as a function of Jh for para- 
magnetic and antiferromagnetic spin configurations for 3 = 
and g = 1.4 . AFM state leads to a reduction in FE dis- 
tortion, b) Phase diagram in A-Jh phase space showing the 
regions of FE stability. 



Via an interplay between charge, spin, lattice and or- 
bital degrees of freedom this leads to the experimentally 
observed magnetically suppressed ferroelectricity. This 
type of strong magneto-electric coupling being present 
in SrLrtBax/aMnOa opens new routes for the search of 
multiferroic materials different from other Mn based ox- 
ides such as RMnC>3 and RMn 2 Os 0], going beyond the 
requirement of the magnetic ordering breaking the inver- 
sion symmetry, whereby it causes a ferroelectric instabil- 
ity. 
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Given a specific configuration of lattice distortions and 
tig spins, one can easily diagonalize the electronic prob- 
lem numerically on finite lattices. We use K s =15 and 
J s =0.1 while calculating the dependencies on the other 
parameters. Our focus is to explore the possibility of a 
FE state in in both the nonmagnetic and AFM phase. 
Therefore, rather than performing a lengthy minimiza- 
tion of the energy as a function of classical spin and lat- 
tice variables we compare the total energy of only the 
relevant configurations. We use the magnitude of off- 
centering distortions u(i) = ito as a variational parame- 
ter and determine the distortions u m in that correspond 
to the lowest total energy. A non-zero value of u m i n is 
the hallmark of a FE state. The onset of antiferromag- 
netism leads to a reduction in the tendency to form a FE 
state, which is reflected in a reduced value of u min for 
AFM order shown in Fig. @]a). 

The results of model calculations are summarized in 
the A—Jh phase diagram in Fig. Hp. The phase diagram 
shows that the FE phase is stabilized over a wide range 
of parameter space when the system is in PM state. The 
presence of AFM order shrinks the regions of stability of 
the FE order, and in general the AFM order reduces 
the value of FE polarization for all parameter values. 
We present the phase diagram for g = 0, which shows 
that the mechanism for FE ordering does not depend on 
the Peierls type electron-lattice coupling. However, the 
presence of a non-zero g further stabilizes the FE order. 

In conclusion, by combining different theoretical ap- 
proaches we highlight the intricate interrelationship be- 
tween magnetic and ferroelectric orderings in rece ntly 
discovered multiferroic phase of Sr^Bai^MnOs [14j | . 
The new mechanism at play relies on the distortion de- 
pendent in-plane hopping between Mn and O sites and 
strongly depends on the onset of the magnetic order. 
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